One-dimensional metallic behavior of the stripe phase in La2_2:Sr2:Cu04 

Marcus Fleck, Alexander I. Lichtenstein and Eva Pavarini 
Max-Planck-Institut fiir Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Federal Republic of Germany 

Andrzej M. Oles 

Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 

(February 1, 2008) 

Using an exact diagonalization method within the dynamical mean-field theory we study stripe 
phases in the two-dimensional Hubbard model. We find a crossover at doping S ~ 0.05 from diagonal 
stripes to vertical site-centered stripes with populated domain walls, stable in a broad range of 
doping, 0.05 < 5 < 0.17. The calculated chemical potential shift oc —S^ and the doping dependence 
of the magnetic incommensurability are in quantitative agreement with the experimental results for 
doped La2-iSriCu04. The electronic structure shows one-dimensional metallic behavior along the 
domain walls, and explains the suppression of spectral weight along the Brillouin zone diagonal. 

PACS numbers; 71.27.-fa, 75.10.-b, 74.72.-h, 79.60.-i. 



It is commonly believed that the understanding of nor- 
mal state properties of high-temperature superconduct- 
ing cuprates (HTSC) will provide important clues for the 
understanding of superconductivity itself. The undopcd 
compounds, La2Cu04 and YBa2Cu306, are insulators 
and exhibit long-range antiferromagnetic (AF) order, 
which is rapidly destroyed and replaced by short-range 
AF correlations as holes are doped into the Cu02 planes 
||,|. Thus, one of the most important features is the 
nature of the interplay between the AF spin fluctuations 
and superconductivity. Incommensurate charge and spin 
order, discovered first in Lai.6-a:Ndo.4Sra;Cu04 sug- 
gests that the strong competition between hole propa- 
gation and AF order in the Cu02 planes leads to segre- 
gation of holes in regions without AF order. These re- 
gions form one-dimensional (ID) substructures, so called 
stripes, which act as domain walls Q. The essentialy 
identical momentum dependence of the magnetic scat- 
tering in La2-a;Sra;Cu04 provides evidence for the 
stripe phases in this class of materials. 

If indeed realized in a broad range of doping, the stripe 
phase should have measurable consequences. Studies of 
La2-a;-yNdj,Sra;Cu04 showed a chemical potential shift 
in underdoped and overdoped cuprates responsible for 
the breakdown of the Fermi liquid picture, a pseudo- 
gap which opens at the Fermi level ||6|, and a real gap 
for charge excitations in the electronic structure around 
momentum (7r/2,7r/2) [0,1). It may be expected that 
such puzzling features follow from strong Coulomb inter- 
actions at Cu ions, and it would be interesting to investi- 
gate whether they are fingerprints of a stripe phase and 
could be reproduced by considering a generic model for 
the HTSC, a two-dimensional (2D) Hubbard model. 

Stripe phases were first found in the Hartree-Fock (HF) g-^-i (k, iu;^) ^ (tuj^ + ^,)5^^ - h,j (k) - Y.,,„ {iuJ^)S^J , (2) 
approximation H, with empty (filled by holes) domain 

walls in an insulating ground state. In contrast, the where /i is the chemical potential, and (k) is an (L x L) 

calculations which include electron correlations indicate matrix which describes the kinetic energy, /iy(k) = 



that the ground state of an AF system with strong short- 
range Coulomb repulsion is a stripe phase with populated 
domain walls at low doping [|0-12|. Hence a partially 
filled band might be expected, and charge transport along 
the walls becomes possible. These results clearly empha- 
size the need for a reliable and controlled approximation 
scheme in order to study the physics of stripe phases. 

In this Letter we present an exact solution of the dy- 
namical mean-field theory (DMFT) [|l3| equations for the 
stripe phase of the 2D Hubbard model. The DMFT ap- 
proach allows to treat the hole correlations in a non- 
perturbative way using local self energy Recently 
we have shown that within DMFT one obtains the cor- 
rect dispersion and spectral weights of quasiparticle (QP) 
states in the Hubbard model at half-filling (n = 1) ||l5|| . 

Here we investigate long-range stripe order in the 2D 
Hubbard model at zero temperature. The square lat- 
tice is thereby covered by N supercells containing L sites 
each, and the ground state energy of the Hamiltonian, 



H = - 



'"ini,n]"'micr"'nj<T 



U 
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has been determined using this constraint. Positions 
Rmi = Tm -|- Vi of nonequivalent sites i — 1, L within 
the unit cell m are labelled by a pair of indices {mi}. 
We focus on the generic behavior of the stripe phase 
and thus restrict the hopping term to nearest-neighbors 
{mi} and {nj} only, tmi.nj — t, and take a uniform on- 
site Coulomb interaction U. The one-particle Green's 
function in the stripe phase is given by an {L x L) ma- 
trix, GijaCk-Tiuj^), on the imaginary energy axis lUi, ~ 
{2v H- l)7rr with fictitious temperature T. It contains a 
site- and spin-dependent local self energy p3|,p^. 



1 



exp(— ik(Roi — 'Rnj))toi,nj ■ The local Green's func- 
tions for each nonequivalent site i are calculated from 
the diagonal elements of the Green's function matrix (|^), 
GiiaiiuJv) = N-^Y.'kGua{\<.,iuj^). Self-consistency of 
site i with its effective medium requires, 



(3) 



similar to the situation in thin films | |l6| |. 

For the solution of the effective impurity model with 
hybridization parameters, Via{k), as well as diagonal en- 
ergies, eia-{k), for each non-equivalent site i in the stripe 
supercell we employed the exact diagonalization method 
of Caffarel and Krauth By fitting Q%^{iuju) on the 
imaginary energy axis. 



k=2 



(4) 



the parameters of an effective DMFT-impurity cluster 
with Us sites are obtained. After solution of the effec- 
tive cluster problem with Lanczos algorithm (ris ~ 8), 
the local Green's function Giia{ioJu), and local electron 
densities were determined. The self-consistency is im- 
plemented by extracting the new selfcncrgy for the next 
DMFT- iteration from Eq. (||). Finally, the Green func- 
tions (H) serve to determine the spectral function. 



A(k,a;) 
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Below we summarize the results obtained for U = 12t, 
a value representative for La2-a;Sra;Cu04 p8[ |, and for a 
broad range of hole doping {6 = 1 — n), 0.03 < 6 < 0.2, 
where we found that the ground state contains populated 
domain walls. First, for 0.03 < S < 0.05, diagonal stripe 
supercells made out of pieces of site-centered vertical do- 
main walls are stabilized by a (weak) charge density wave 
(CDW) superimposed with a spin density wave (SDW) 
along the wall. The SDW domain wall unit cells con- 
sist of four sites, |0) - | T) - |0) - | i). For 6 = 0.05 
we found the somewhat lower electron densities at non- 
magnetic (|0)) sites {uq ~ 0.848) than the densities at 
magnetic (|cr)) sites (rim — 0.860) with magnetic mo- 
ment of TO ~ 0.334. The SDW unit cell is stabilized by 
only ~ 0.584 holes, which indicates that such states are 
precursors of the undoped AF Mott insulator [p^ . 

In agreement with neutron scattering experiments , 
we found a particular stability of the vertical stripes with 
populated domain walls. As the most robust structure of 
the doped systems with 0.05 < S < 0.17 we identified the 
site-centered stripe phase (Fig. pj). This phase is stabi- 



lized by electron correlations |10|-|12| and by the kinetic 
energy gains on the populated domain walls which desta- 
bilize a SDW along the walls in this doping regime. At 
higher doping S > 0.17 the bond-centered stripe phase of 
White and Scalapino mU is energetically favored, kinks 



and antikinks along the domain walls occur, and the 
stripe structure gradually melts. 
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FIG. 1. Vertical stripe phases for 5 = 1/12 (left) and 
S — 0.15 (right) at i7 = 12t. Top part shows doped hole 
(circles) and magnetization density (arrows). The spatial 
variation of these quantities are represented by Uhilx) (open 
squares) and STr(lx) (solid circles) in the lower part. 

The size of AF domains in the site-centered stripe 
phase shrinks with increasing doping for (5 < 1/8 (Fig. 
|l|). For example, the charge unit cell contains eight 
(four) sites at doping 5 = 1/16 {5 = 1/8), while the 
electron density on the sites of the walls is almost con- 
stant, being ni ~ 0.850 and rii ~ 0.830, respectively. 
In agreement with the results of the numerical den- 
sity matrix renormalization group pT| , the self-consistent 
DMFT densities in the stripe unit cell are character- 
ized by much smoother variations than in the corre- 
sponding HF states Beyond 5 — 1/8 we find a 
lock-in effect of the same structure with a charge (mag- 
netic) unit cell consisting of four (eight) sites, and the 
doped hole density, n/,(l^) = 1 - (n-(i,,o),T + "-(i=.,o)a)> 
increasing faster within the AF domains than on the 
wall sites (Fig. The magnetic domain structure is 
best described by the modulated magnetization density, 

SAlx) = ij;^Ei„(-i)'"=+'4('^(;x^).T ^ "^(i.Ma)^ pro- 
jected on the direction perpendicular to the wall [ pT| . 

The stability of the above stripe phases is investigated 
by the ground state energy normalized per density of 
doped holes. Eh = [Eo{5) - E{Q)]/5, where Eo{5) is the 
ground state energy at doping 5. The energy Eh/t is a 
monotonically increasing function of doping (Fig. ||(a)), 
showing that the different stripe phases are stable against 
macroscopic phase separation. From our results we con- 
clude that short-range Coulomb repulsion suffices to ob- 
tain populated domain walls over a wide range of dop- 
ing. The diagonal stripes, stable at low doping 5 < 0.06, 
are followed by vertical site-centered stripes which com- 
pete with bond-centered stripes, both having a consid- 
erably larger energy gain than homogenous phases, such 
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as spin spirals. We have verified that the kinetic energy 
is gained mainly on the sites which belong to the do- 
main walls. Such energy gains are larger at the nonmag- 
netic domain walls (Fig. ||(a)) than on the magnetic sites 
of bond-centered domain walls. At the same time the 
small energy difference between these both quite differ- 
ent states indicates a strong tendency towards transverse 
stripe fluctuations which might enhance superconducting 
correlations in the ground state [ p^ . 

We have found that the chemical potential shifts down- 
wards with hole doping, oc —5"^ (Fig. ||(b)), in agree- 
ment with the experimental results of Ino et al. Q , and 
with the Monte-Carlo simulation of a 2D Hubbard model 
pO| . Therefore, the charge susceptibility is enhanced to- 
wards 5 —t 0, reproducing a universal property of the 
Mott-Hubbard metal- insulator transition |20|. 
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FIG. 2. Evolution of the stripe phases with doping 5 
{U = 12t): (a) Energy per doped hole Eh/t in diagonal 
SDW stripes (squares) , and vertical (circles) site-centered and 
bond-centered (triangles) stripes; filled (empty) symbols show 
the energy in the ground (excited) states, (b) Shift of the 
chemical potential A/x/t (filled symbols) and the quadratic 
fit Af/,/t = aS'^, with a = -15.57 (dashed line), (c) Shift 
ri of the maxima of ^(Q) in: diagonal SDW stripes with 
Q = [(1 ± \/2T])'rr, (1 ± y/2ri)Tr] (filled squares), and vertical 
site-centered (filled circles) and bond-centered (filled trian- 
gles) stripes, both with Q = [(1 ± 2ri)Tv,TT]. Empty symbols 
show the da ta p oints of Yamada et al. M (circles) and Waki- 
moto et al. pil (squares). 

As both the bond-centered and site-centered stripe 
phases have the same size of the magnetic unit cell, they 
give the same pattern in neutron scattering and are thus 
indistinguishable experimentally. The neutron scatter- 
ing structure factor S'(Q) in the stripe phase has the 
maxima shifted away from the M = (tt, tt) point to 
Q = [(1 ± 2?7vert)7r, tt] points for the structures of Fig. 
^ [and to Q = [tt, (1 ± 277vGrt)'''] for equivalent horizontal 
stripes]. The present calculations give a linear depen- 
dence revert = S for S < 1/8 and ryvort = 1/8 for 5 > 1/8 
(Fig. 0(c)). Such a behavior was observed by Yamada et 



al. Q , and indicates a unique stability of populated do- 
main walls in the stripe phase. The correlations included 
within the DMFT and its capability to describe the Mott- 
Hubbard metal-insulator transition play thereby a crucial 
role, as other filling and periodicity of the stripe phase 
are found in HF calculations Also the points found 
at low doping (Fig. ||(c)) corresponding to the diago- 
nal stripe structures Q = [(1 ± 2?7diag)7r, (1 ± 2?7diag)7r] 
agree perfectly well with the recent neutron experiments 
of Wakimoto et al. We find ''7diag — 5/ %/2, where 

the factor 1/^/2 is due to the rhombic lattice constant 
in diagonal stripe structures as suggested by experiment 
|2l| . This results in the relation f^diag = ??vcrt /%/2, and 
the linear dependence revert — 5 holds also for 5 < 0.06 
(Fig. H(c)). Furthermore, our calculations predict verti- 
cal SDW domain wall unit cells in the diagonal stripe 
phase. Thus additional elastic magnetic superlattice 
peaks should be visible in neutron scattering experiments 
around Q = [(1 -I- 27ydiag)7r/2, (H- 277diag)7r] with a weight 
smaller by a factor ~ 3.7, if such phases do exist in heav- 
ily underdoped La2-a;-yNdySra;Cu04. 
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FIG. 3. Spectral function A(k, oj) along the main directions 
of the 2D BZ of the stripe phase at 5 = 1/12 with U = 12t. 

Let us focus on the spectral function A(k, lS) of the 
stripe phases shown in Fig. |l|. The photoemission w < /i 
spectra consist of the lower Hubbard band at an energy 
Lo — [i — 4.8i and low-energy states, well separated from 
the Hubbard band, extending over an energy range of 
^ 2t. The QP states known from the dispersion of a 
single hole in the t-3 model survive also in the stripe 
phase up to (5 = 0.15 and are characterized by a con- 
siderable spectral weight and a bandwidth ~ 2J (here 
J/t = /It/U = 1/3) [H (Figs. I and I). Due to the stripe 
structure we find that the directions Y — X [X = (tt, 0)] 
and Y — Y [y = (0, tt)] are nonequivalent. 

At (5 = 1/12 we observe a pseudogap at the X point 
(Fig. |). The QP weight there is composed out of QP 
states originating from the dressing of a moving hole by 
quantum fluctuations in an AF background and local- 
ized states from the ID electronic structure of the site- 
centered stripe phase This superposition of QP 
weight explains the flat hand around the X point and the 
Fermi level crossing at (tt, 7r/4) observed in recent angle- 
resolved photoemission (ARPES) experiments On 
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the contrary, the QP's originating from the ID features 
of the stripe phase, are not seen in ARPES around the Y 
point, as the structure factor vanishes [^, and one only 
resolves the spin-polaron QP with dispersion ^ 2J. 

As the most spectacular result, a gap for charge ex- 
citations opens in the underdoped regime at the Fermi 
energy around S = (7r/2,7r/2) (Fig. ^. Little or no 
spectral weight is found at momenta k = (7r/4,7r/4) and 
k — (0,7r/4) vsfhere, notably, the ID electronic structure 
should show Fermi level crossings |^^. This behavior 
agrees quantitatively with the ARPES measurements on 
La2-:,Sr^Cu04 and Lai.28Ndo.6Sro.i2Cu04 |]. 

In order to understand the gap structure we calculated 
the electronic structure of H = —tJ2<ij> a ^i<7^j<j + 
"^iUiSf, where Ui are site-dependent spin potentials in 
the stripe supercell. Local energy contributions Vi = 
Ui\{Sf)\ are treated as parameters, where Vb = on 
the domain wall. Wc find a condition for vanishing 
photoemission weight at momentum k = (7r/4,7r/4), 
Vo-(-2/^ > 2Vb+/^, which is satisfied by the present self- 
consistent DMFT spin densities, with — 2.07i and 
Vb+;^ — 0.99t (Fig. The strong renormalization of Vi 
in DMFT is due to charge fluctuations and demonstrates 
that local correlations play a crucial role in understand- 
ing the ARPES spectra of HTSC gT 
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FIG. 4. A(k,a;) as in Fig. |but at 5 = 0.15. 

The ARPES spectral weight at (5 = 1/12 around 
uj — fi ^ —2t entirely originates from bands of the stripe 
supercell (Fig. ||). As expected, spectral weight is trans- 
ferred with increasing doping from that energy region to 
the inverse photoemission w > /x (Fig. ^) . Finally, we ob- 
served that the gaps at the Y and S point are gradually 
filled by spectral weight as the stripe order melts. 

Summarizing, we have shown that vertical stripe 
phases with populated domain walls are robust structures 
in a broad range of S. Their spectral properties show an 
interesting superposition of the QP's known from doped 
2D antiferromagnets with a ID metallic behavior. Such 
experimental features as: (i) the chemical potential shift 
A/x oc —S^, (ii) the incommensurability of spin fluctu- 
ations, and (iii) the gradual disappearence of the pho- 
toemission flat band and the pseudogap (gap) at the X 
(S) point, find a natural explanation and accompany a 
gradual crossover from the stripe phases into a (strongly 



correlated) Fermi liquid with increasing doping. 
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